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Hadron masses are calculated in quenched lattice QCD on a variety of lattices in order to probe the 
scaling behavior of the Fat-Link Irrelevant Clover (FLIC) fermion action, a fat-link clover fermion 
action in which the purely irrelevant operators of the fermion action are constructed using APE- 
smeared links. The scaling analysis indicates FLIC fermions provide a new form of nonperturbative 
0(a) improvement where near-continuum results are obtained at finite lattice spacing. 
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Discretization of the continuum action of QCD for use 
on a space-time lattice grid with finite lattice spacing 
a introduces errors of order a™. For Wilson fermions 
n = 1 due to the introduction of an irrelevant energy 
dimension-five lattice-Laplacian operator designed to re- 
move fermion doublers from the naive lattice theory by 
giving the doublers a mass proportional to a -1 . While 
Wilson fermions are computationally inexpensive, the ap- 
proach to the continuum limit is slow. Observables are 
spoiled by large 0(a) discretization errors and quantita- 
tive extrapolations to the continuum limit must be per- 
formed from simulations performed at small lattice spac- 
ings typically less than 0.1 fm. This necessitates the use 
of very large four-dimensional lattices in order to pro- 
vide reasonable physical simulation volumes exceeding 2 
fm on a side. 

Since the cost of simulations increases as a~ n with n 
exceeding 5, considerable savings have been achieved by 
designing lattice fermion actions free of 0(a) discretiza- 
tion errors. A systematic approach to achieving Oia) 
improvement of the lattice fermion action in general |3] is 
to consider all possible gauge invariant, local dimension- 
five operators, respecting the symmetries of QCD 
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Operator 0\ is a new local operator in the lattice fermion 
action and must be included. On the other hand, 03 and 
05 of Eq. Q act to simply renormalize the coefficients 
of existing terms in the lattice action, removing 0(am q ) 
terms from the relation between bare and renormalized 
quantities 0]. 

The key observation to efficient 0(a) improvement is 
that the 0(a) improvement afforded by two- link terms of 
the fermion action may be incorporated to 0(a) into 
the standard Wilson fermion action complemented by 0\ 



though the following transformation of the fermion fields 

ip — * ip' = (1 + b q r am q ) (1 — c q r a B) ip , 

ip ^ ip' = (l + b q ram g )ip(l + c q ra^) , (2) 

where ip' represents the physical fermion field recovered 
in the continuum limit, while tp is the lattice fermion 
field used in the numerical simulations. At tree-level, 
b q = c q = 1/4 correctly incorporates the 0(a) corrections 
of 02 and 04 into the fermion action. 

In summary, 0i, the "clover" term, is the only 
dimension-five operator explicitly required to comple- 
ment the Wilson action to obtain 0(a) improvement. 
This particular action is known as the Sheikholcslami- 
Wohlert fermion Q action 
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where Sw is the standard Wilson action jfj, and Csw is 
the clover coefficient which can be tuned to remove 0(a) 
artifacts to all orders in the gauge coupling constant g. 

While this action has been known for some time, the 
difficulty has been in accurately determining the renor- 
malization of the improvement coefficients, Csw, b q , c q , 
etc., in the interacting quantum field theory. At the lat- 
tice spacings typically considered in today's lattice sim- 
ulations, the tree level values can differ by a factor of 
two from the renormalized values. While mean-field im- 
proved estimates of the renormalized coefficients provide 
substantial corrections, they are not sufficiently accurate 
to remove 0(a) errors to all orders in the gauge coupling 
constant g 0- Nonperturbative (NP) 0(a) improvement 
[U tunes Csw to all powers in g 2 and displays excellent 
scaling as discussed further in the following. 

In a previous paper 0, we introduced a new form of 
0(a) fermion-action improvement in which the renormal- 
ization of improvement coefficients is addressed in a very 
different manner. Central to the approach is the obser- 
vation that the fermion doublers of the naive theory are 
removed by the Wilson term at tree level. In place of ap- 
plying techniques to estimate the renormalization of the 
improvement coefficients induced by the gauge fields of 
QCD, techniques are applied to modify the gauge fields 
to suppress renormalizations such that tree-level knowl- 
edge of improvement coefficients is adequate. 
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There are many accepted methods for removing short- 
distance fluctuations from gauge field configurations in- 
cluding APE smearing, HYP-smearing, and their vari- 
ants. The central feature of FLIC fermions is to construct 
the fermion action using two sets of gauge fields. In the 
lattice operators providing the relevant dimension-four 
operators of the continuum action, one works with the 
untouched gauge fields generated via Monte Carlo meth- 
ods, while the smoothed gauge fields are introduced only 
in the purely irrelevant lattice operators having dimen- 
sion five or more. We refer to this action as the Fat-Link 
Irrelevant Clover (FLIC) fermion action. 

The motivation behind constructing the FLIC action 
is to benefit from the reduced exceptional configuration 
problem of fat-link actions while retaining short- 
distance quark interactions in the relevant operators of 
the fermion action. The expectation is that improvement 
in the condition number of the FLIC fermion matrix will 
allow rapid calculations of fermion propagators and effi- 
cient access to the chiral limit of full QCD ^(| • 

In this letter we present the first comprehensive scal- 
ing analysis of FLIC fermions where four different lattice 
spacings are considered on five lattices each sampled by 
200 configurations. The scaling analysis shows convinc- 
ingly that the new FLIC fermion action removes 0(a) 
lattice artifacts from the Wilson fermion action, on a 
level comparable to that of the nonperturbative improved 
clover action. 

The simulations are performed using an C(a 2 )-mean- 
field improved Luscher-Weisz plaquette plus rectangle 
gauge action [U on 12 3 x 24, 16 3 x 32 and 20 3 x 40 lattices 
with lattice spacings of 0.093, 0.122, 0.134 and 0.165 fm 
determined from a string tension analysis incorporating 
the lattice coulomb potential ^3 with ^fa = 440 MeV. 
Initial studies of FLIC, mean-field improved clover and 
Wilson quark actions were made using 50 configurations. 
The scaling analysis of FLIC fermions presented here 
is performed with a total of 200 configurations at each 
lattice spacing and volume. Gauge configurations are 
generated using the Cabibbo-Marinari pseudo heat-bath 
algorithm with three diagonal SU(2) subgroups looped 
over twice. Simulations are performed using a parallel 
algorithm with appropriate link partitioning |l3j . and 
the error analysis is performed by a third-order, single- 
elimination jackknife, with the % 2 per degree of freedom 
(Wdf) obtained via covariance matrix fits. 

Fat links 0, 0] are created by averaging or smearing 
links on the lattice with their nearest transverse neigh- 
bors in a gauge covariant manner (APE smearing). The 
smearing procedure 0] replaces a link, U^x), with a 
sum of the link and a times its staples 



U^(x)^U' fl (x) = (l-a)U fl (x) (4) 
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followed by projection back to SU(3). We select the uni- 
tary matrix C/J L which maximizes 



■Retv(U^U'J), 



(5) 



by iterating over the three diagonal SU(2) subgroups of 
SU(3). Performing eight iterations over these subgroups 
gives gauge invariance up to seven significant figures. 
The combined procedure of smearing and projection is 
repeated to create a fat link. 

The mean-field improved FLIC action now becomes 
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where is constructed using fat links, Uq is the mean 
link calculated with fat links, and where the mean-field 
improved Fat-Link Irrelevant Wilson action is 
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with k — l/(2m + 8r). We take the standard value r = 1. 
Our notation uses the Pauli (Sakurai) representation of 
the Dirac 7-matrices defined in Appendix B of Saku- 
rai [20|. In particular, the 7-matrices are hermitian and 
cr^ = [7 M , 7^]/(2«). 

As reported in Ref. [l|, the mean- field improvement 
parameter for the fat links is very close to 1. Hence, the 
mean-field improved coefficient for Csw is expected to 
be accurate. A significant advantage of the fat-link irrel- 
evant operator approach is that one can now use highly 
improved definitions of F^ v which give impressive near- 
integer results for the topological charge |21j . 

In particular, we employ the 3-loop O(a 4 )-improved 
definition of F^ v in which the standard clover-sum of four 
lxl loops lying in the /i, v plane is combined with 2x2 
and 3x3 loop clovers. Bilson-Thompson et al. [2lJ find 
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where W nxn is the clover-sum of four n x n loops and 
F^v is made traceless by subtracting 1 /3 of the trace from 
each diagonal element of the 3x3 color matrix. This defi- 
nition reproduces the continuum limit with 0(a 6 ) errors. 
On approximately self-dual configurations, this operator 
produces integer topological charge to better than 4 parts 
in 10 4 . We have also considered a 5-loop improved F^ 
which agrees with the 3-loop version to better than 4 
parts in 10 4 HJj. 
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FIG. 1: Nucleon and vector meson masses for the Wil- 
son, Mean-Field (MF) improved, NP-improved clover, domain 
wall, fixed point, improved staggered and FLIC actions ob- 
tained by interpolating simulation results to m^/m p = 0.7. 
Results from the current simulations are indicated by the solid 
symbols; those from earlier simulations by open or hatched 
symbols. The solid-lines illustrate fits, constrained to have a 
common continuum limit, to FLIC, NP-improved clover and 
Wilson fermion action results obtained on physically large lat- 
tice volumes. Further details are provided in the text. 

An important consideration is the amount of smearing 
to apply to the gauge fields of the irrelevant operators. 
Since the aim is to remove perturbative renormalizations 
of the improvement coefficients, we monitor the 3-loop 
0(a )-improved topological charge constructed with the 
3-loop 0(a 4 )-improved definition of F^ v of Eq. JSJl as 
a function of smearing sweep. The topological charge 
is known to have a large multiplicative renormalization 
|19| and serves as an ideal operator for monitoring the 
removal of perturbative physics under smearing. We find 
that the topological charge varies rapidly over the first 
few sweeps of smearing but then makes only small varia- 
tions thereafter. We define the optimal number of sweeps 
to be the minimum number of sweeps required to reach 
the smoothly varying regime. It is interesting that our 
findings for optimal smearing coincide with that required 
to provide the optimal condition number of the FLIC 
fermion matrix in the negative mass regime relevant to 
overlap fermions |17| . 

Hadron masses are extracted from the Euclidean time 
dependence of the calculated two-point correlation func- 
tions using standard techniques [24|. To compare with 
the results of Ref. Q, we interpolate our results to a 
pseudo-scalar to vector meson mass ratio of m ff /m p = 
0.7. The scaling behavior of the various fermion actions 
is illustrated in Fig. ^ Lattice volumes, spacings and 
hadron masses in units of the string tension are given in 



TABLE I: Lattice parameters and results for the vector meson 
and nucleon masses interpolated to m^/mp = 0.7. 
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m v l\fa 




MO 


4.38 


16 3 


x 32 


200 


0.371 


2.360(20) 


3.439(27) 


0.8761 


4.53 


20 3 


x 40 


200 


0.299 


2.324(15) 


3.382(24) 


0.8859 


4.60 


12 3 


x 24 


200 


0.274 


2.434(26) 


3.554(33) 


0.8889 


4.60 


16 3 


x 32 


200 


0.274 


2.336(22) 


3.400(26) 


0.8889 


4.80 


h. 


x 32 


200 


0.210 


2.427(23) 


3.538(61) 


0.8966 



Table H 

Since the smearing radius |16| is proportional to the 
product of a and the number of smearing sweeps, n, we 
fix a = 0.7 and vary n. For the fine lattices with a 2 a ~ 
0.075 and 0.045, four smearing sweeps are performed. For 
the coarser lattices with a 2 a ~ 0.09 and 0.14 we perform 
six and eight smearing sweeps respectively. The effective 
range of the smearing [Til within which interactions are 
suppressed is (r 2 ) = a 2 an/3, providing an RMS radius 
of 1.0, 1.2 and 1.4 lattice spacings for n = 4, 6 and 8 
smearing sweeps at a — 0.7 respectively. 

Actions with fat-link irrelevant operators perform ex- 
tremely well, lying very near the horizontal dashed lines 
corresponding to continuum limit results at finite lattice 
spacing. For reference we have also calculated masses for 
the Wilson action at a 2 a ~ 0.075 which agree with those 
of Ref. [7j . We also compare our results with the standard 
Mean-Field Improved Clover (MFIC) action. We mean- 
field improve as defined in Eqs. © and (JJJ) but with thin 
links throughout. For this action, the standard 1-loop 
definition of is used. For both the vector meson and 
the nucleon, the FLIC actions perform significantly bet- 
ter than the mean-field improved clover action. 

Finally, our FLIC results compare extremely well with 
a variety of improved actions found in the literature. In 
particular, in Fig. ^ we compare with results using do- 
main wall |25l. fixed point |26j| and improved staggered 
fermions [27]]. 

The two different volumes used at a 2 a ~ 0.075 reveal 
a small finite volume effect, which increases the mass for 
the smaller volumes at a 2 a ~ 0.075 and ~ 0.045. Exam- 
ination of points from the small and large volumes sepa- 
rately indicates continued scaling toward the continuum 
limit. While the finite volume effect will produce a dif- 
ferent continuum limit value, the slope of the points from 
the smaller and larger volumes agree, consistent with er- 
rors of 0{a 2 ). 

Focusing on simulation results from physical volumes 
with extents ~ 2 fm and larger, we perform a simulta- 
neous fit of the FLIC, NP-improved clover and Wilson 
fermion action results. The fits are constrained to have 
a common continuum limit and assume errors are 0{a 2 ) 
for FLIC and NP-improved clover actions and 0(a) for 
the Wilson action. To obtain a data point at our fine 
lattice spacing, a 2 a ~ 0.045, we use the observed fi- 
nite volume effect at a 2 a ~ 0.075 to correct the point 
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TABLE II: Fit parameters and xof f° r joint and separate fits 
of the FLIC, NP-improved clover and Wilson hadron masses 
We fit to an ansatze of the form mn Isftj = Ho + Hi a^/a + 
H2 a 2 a, where the hadron, H, can be the vector meson, V , or 
the nucleon, N. 





FLIC 


NP Clover 


Wilson 






v 2 


v 2 


Vi 


v 2 


2 

Xdf 


2.324(24) 


0.18(23) 


0.69(18) 


-1.78(20) 


0.17(43) 


0.96 


2.317(18) 


0.24(18) 


0.74(14) 


-1.71(7) 





0.87 


2.320(25) 


0.22(24) 


0.72(19) 






0.76 


2.291(34) 


0.47(32) 








0.40 


No 


N 2 


N 2 


Ni 


iV 2 


2 

Xdf 


3.415(55) 


-0.05(56) 


0.81(41) 


-2.18(45) 


0.34(97) 


1.88 


3.402(39) 


0.07(42) 


0.90(31) 


-2.04(15) 





1.69 


3.402(50) 


0.08(51) 


0.90(38) 






1.52 


3.335(53) 


0.71(53) 








0.39 



at a 2 a ~ 0.045 as illustrated by the open triangles. The 
results from these fits are given in Table ITT1 In light of 
the different gauge actions used in the analyses and the 
fact that the lattice volumes considered are not perfectly 
matched throughout all the simulations, an acceptable \ 2 
per degree of freedom is obtained for both the nucleon 
and p- meson fits. 



To assess the sensitivity of our results on the num- 
ber of smearing sweeps, we perform a second calculation 
with only four smearing sweeps for our lattice having the 
coarsest lattice spacing. This result is indicated by the 
solid diamond symbol offset to the right for clarity. These 
two results at a?o ~ 0.135 reveal an insensitivity to the 
number of APE-smearing sweeps used in constructing the 
irrelevant operators of the FLIC fermion action. This in- 
sensitivity suggests that one could define the FLIC ac- 
tion in terms of a fixed number of APE smearing sweeps 
independent of the lattice spacing. Upon taking the con- 
tinuum limit the smearing radius would still tend to zero 
as required. 

In conclusion, the use of fat links in the irrelevant op- 
erators of the FLIC fermion action provides a new form 
of nonperturbative 0{a) improvement without the need 
of nonperturbative fine tuning. In addition, the 0(a 2 ) 
errors are small for this action. FLIC fermions display 
nearly perfect scaling, providing near continuum limit re- 
sults at finite lattice spacing. 
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